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ABSTRACT 


Experimental data are combined with one-dimensional 
conservation relations to yield information on the energy 
deposition ratio in a parallel-plate accelerator, where 
the downstream flow is confined to a constant area channel. 
Approximately 70% of the total input power is detected in 
the exhaust flow, of which only about 20% appears as 
directed kinetic energy, thus implying that a downstream 
expansion to convert chamber enthalpy into kinetic energy 
must be an important aspect of conventional high power 
MPD arcs. 

Spectroscopic experiments on a quasi-steady MPD argon 
accelerator verify the presence of AIII and the absence of 
AI , and indicate an azimuthal structure in the jet related 
to the mass injection locations. Measurements of pressure 
in the arc chamber and impact pressure in the exhaust jet 
using a piezocrystal backed by a Plexiglas rod are in good 
agreement with the electromagnetic thrust model. 

Other experiments related to various aspects of MPD 
accelerator operation continue: additional characteristics 
of a hollow-cathode MPD arc are determined; first indica- 
tions of thrust-time profiles are obtained; a study is 
initiated on the biased double probe in the exhaust plasma; 
and mass pulses are recorded from a new turbo- in j ec tor 
device . 
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I . INTRODUCTION 


This semi-annual report contains a review of a 
Ph.D. research program recently completed in our labora- 
tory, along with more concise descriptions of several 
other programs presently in progress. The first section 
outlines the use of a variety of experimental results in 
a one-dimensional model to determine the energy deposi- 
tion ratio in our parallel-plate accelerator, a major 
portion of the Ph.D. thesis of M. S. DiCapua. In the 
second section, results from various studies on the 
pulsed MPD accelerator are presented, including velocity 
and species data based on both time-resolved and time- 
integrated spectroscopic techniques and the use of a 
piezoelectric pressure probe to determine the chamber 
pressure and impact pressure profile in the exhaust. 

The final section reviews the progress of other programs 
in various stages of completion from preliminary and 
feasibility tests to well established areas of research. 
These include terminal data on an MPD arc with a hollow 
cathode whose outer surface is insulated, further progress 
with the time-resolved thrust measurement technique, a 
study to determine the applicability of the biased double- 
probe as a diagnostic tool in the MPD exhaust, and the 
first tailored mass pulses generated with the turbo-pulse 
gas injector. More detailed reports on these programs, as 
well as on other studies in progress will be deferred 
until the following report. 
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II. ENERGY DEPOSITION IN THE PARALLEL- PLATE 
PLASMA ACCELERATOR <DiCapua) 

A. Introduction 

The central interest of this study is the initial appor- 
tionment of the input power to the streaming and thermal 
energy of the flow in a plasma accelerator. This initial 
energy deposition ratio achieves great importance if one is 
committed to the task of efficient plasma acceleration. In- 
sufficient detailed knowledge of the electromagnetic, gas- 
dynamic, and thermodynamic fields throughout the acceleration 
region commits us to obtain the energy deposition ratio 
through a black box or terminal analysis of the accelerator. 
This analysis is based upon measurement of the terminal electr 
cal parameters of the accelerator and upon observation of prop 
erties of the flow into and out of the acceleration region. A 
self-consistent picture can then be drawn through the applica- 
tion of momentum and energy conservation arguments. Even on 
this basis, however, several complications arise due to the 
complex non-equilibrium state of the working fluid and the in- 
troduction of spurious mass into the flow through electrode 
and insulator erosion as well as ingestion of mass associated 
with the gas.dynamic flow field established prior to the ini- 
tiation of the discharge. 

The parallel plate geometry is particularly well suited 
for this analysis. The electromagnetic fields upstream of the 
acceleration zone can be obtained through simple one dimen- 
sional arguments. Downstream of the acceleration zone, the 
plasma flow field is amenable to a steady, one-dimensional 
treatment as well. The first task, therefore, is to obtain 
the conservation relations .in a convenient form for the 
analysis . 

B. The Conservation Relations 

Analysis of the problem in a terminal fashion can best 
be performed with the conservation relations cast in integral 
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form. This approach allows treatment on the basis of mass, 
momentum, and energy fluxes across appropriate surfaces 
without requiring any detailed knowledge of the processes 
which occur within the acceleration region. 

1. Mass Conservation 

The mass conservation equation for an arbitrary fixed 
volume can be written as: 11 ^ 

K Jv 

where f is the density, u , the velocity, "V, the total 
volume# A, the total surface area of that volume, and dV and 
dA are respectively the volume and surface elements. 




Momentum Conservation 

The momentum conservation equation in integral form can 
II-l 


be written as : 


JA ( 2 ) 


j'L(fu)Jv + r fu(u-jx) _ -jpJA 4 Jfjv 

V A A V /\ 

where p is the pressure, f , the body force per unit volume, 

and Z / viscous stress tensor. In this case, where elec- 

-k. 

tromagnetic forces are present, f takes the form: 

/ * f X 3 (3) 

where j is the current density and B is the magnetic 
field. 

To transform the volume integral of the body forces into 

the required surface integral, the body force must be cast in 

II-2"' 

the form of the divergence of a tensor. The tensor 

can then be integrated over the surface. To this end, the 
body force becomes: 

f _» £ * 3 _= ~ (v *8 ) * 3 (4) 


through the use of Ampere's Law in the low frequency limit: 

(5) 


1 7x3 = 


/•/ 
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Using the vector identity: 

(Vx S) x S = (s'- v)&-{ v(S-B) 

with: . 

o =■ & (v-&) 

added to both sides, the body force becomes; 

fx B = 2 (Vxs)x a = 7-x /as -4 

" r° ^ o 2 


( 6 ) 


(7) 


( 8 ) 


where 


3=3-3 

Z - in 


( 9 ) 


And J- - 6 lj 

The right hand term of equation (8) is the divergence of 
that part of the Maxwell stress tensor which is associated 
with the magnetic field. 

By Gauss' s theorem the volume integral is transformed 
into the desired surface integral: 

I f 1 

( 10 ) 

7 ' 

4 


/ fjv = ^ / 6 (b-ja) -m j 3 ‘ja 


V' A 

and the momentum equation becomes : 


/ Jt + ffu(u-JA)-- -$(p+£)jA+£B(BJA)+ft.jA- O') 

A a ' /A a 


V A 

3. Energy Conservation 


The energy conservation equation in integral form is: 

J (t '+ -)] JV +J « - / pu-JA 

{(u-v).Ja £ f fe)jv f ff.JX 


II-l 


+ 

A 

where 


( 12 ) 


A ) 3* 

is the internal energy per unit mass, rr 


is the 


it */*• 

time rate of change of energy per unit volume as sociated .with 
the magnetic field, and q is the heat flux vector. The 
quantity j. E where j is the current density as before 
and E is the electric field, represents the rate of con- 
version of electromagnetic field energy into mechanical and 
thermal modes. 
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It is again advantageous to transform the volume inte- 
grals associated with the electromagnetic fields into surface 
integrals which are more closely related to the terminal 
parameters of the accelerator. To this end, using Ampere's 
Law: 

B* l 7 c — d E>‘ 


a & + r e = — a 2 

dt 2^' ° dt z/o 


Z2L§.? 


With the identity: 

- e - (v*&) = v- (s’* b) 


and Faraday's Law: 

VX? 

equation (13) becomes: 

A &z 

dt fo 

Thus, using Gauss's theorem. 


= 

Jt 


A ~ - V- (e* b) 

At z yj c " ' ' 


(13) 


(14) 


(15) 


(16) 


f (lir,. - / 


a*B) 


/*• 


Ja 


(17) 


is the Poynting vector which represents the 


The term — — 

flux per unit area per unit time of electromagnetic energy. 
The energy equation therefore takes the form: 


v 


+ f (“■?)• Ja + / T-Ja 

A A 7 0 

where the enthalpy fa is defined as: 


(18) 


fi - e + f/f 

C. The Steady One-Dimensional Approach 


(19) 


The conservation relations derived- in- the -previous sec- - 

tion can be greatly simplified under the assumption of steady 
one-dimensional flow. With this assumption all terms contain- 
ing time derivatives vanish and all properties of the flow 
and the electromagnetic field become the function of only one 
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coordinate, namely the x coordinate as shown in Fig. II-l. 
The relations are further simplified with the assumption of 
negligible viscous and heat transfer effects to be justified 
later. 

The geometry of the accelerator is shown schematically 
in Fig.H-1. The channel has an area A of width w and 
height h. The injector orifices in the back wall have a 
total area A^. The injected argon propellant has a pres- 
sure at the exit plane of the orifices of , a density , 

and a velocity u^ in the x direction. A current J flows 
through the backstrap and accelerator causing a voltage drop 
across the electrodes of V. Downstream of the acceleration 
region the pressure is p,the density f , and the velocity u 
in the x direction. 


Under the steady, one-dimensional assumption, the elec- 
tromagnetic fields are obtained from the current and voltage 
across the accelerator. Although the voltage drop across the 
terminals is \J, only a portion of this is imposed upon the 
plasma flow due to the anode and cathode fall voltages. An 
empirical correction factor oC is used to denote that fraction, 
i.e oC V J represents the fraction of the total input power 
available to the flow. It follows. that the electric field 
everywhere within the plasma except in the falls can be de- 
termined from 

Jf t -J£ = (20 > 

from which 

£ = £ ( 21 ) 

assuming E-j_ , the electric field as defined above, is con- 
stant along the height. By Ampere's Law, on the other hand, 

/ &■ Jl = A J (22) 

and therefore the magnetic field is: 






J 


2 . 


( 23 ) 


(V 



DISCHARGE REGION 
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FIGURE H-l 

AP 25- 4730 


SCHEMATIC OF PARALLEL PLATE ACCELERATOR 
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The Conservation relations therefore become: 

f, A, — f Li A Mass (24) 

fiSA - r, <A, = f>A, + A (pc - f + Ajf) Momentum (25) 

fu A (ft + ^ J ~ f t U' A, y±J sz oC\/J Energy (26) 

D. Experimental Results 

The experimental data strongly suggest that the steady 
analysis is applicable to this accelerator. A typical simul- 
taneous record of current J, terminal voltage V and in- 
ductive voltage Vq taken at the 85 kA, 96g/sec operating 
condition is shown in Fig. II-2a.The current rises in 10 jisec 
and remains constant for 180 jisec. The inductive voltage, 
which is associated with changes of magnetic flux in the 
accelerator due to time derivatives of the total current or 
motion of the current distribution, drops to a negligible 
value after 50 /isec and remains practically zero until 180 jisec 
when the current begins to drop. The voltage V, which is 
measured downstream of the current distribution, represents 
the sum of the resistive voltage and the induced emf across 
the accelerator. This voltage assumes a constant value after 
50 ^tsec and begins to drop at 180 ^isec, indicating steady 
electrodynamic operation of the accelerator ' in this time in- 
terval . 

Measurements of pressure in the accelerator channel show 
that steady gasdynamic operation is achieved as well. Fig. II- 2b 
shows triple overlays of simultaneous pressure measurements 
taken 10 cm upstream of the electrode insulator junction and 
23 cm downstream of the junction in the centerline of the 
channel top. The peak in the upstream pressure p c is 
associated with' the unsteady phase of the discharge. The 
second and third peaks are due to internal reflections in 
the probe so that the upstream pressure achieves a low but 
steady value of about 0.03 atm roughly at the same time the 
voltage reaches the steady phase. 
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The downstream pressure, on the other hand, achieves a 
steady value at 95 jisec and remains roughly steady at approx- 
imately 0.45 atm until 270 psec. There is a time delay be- 
tween achievement of the steady phase at the electrode in- 
sulator junction and at the downstream station due to the 
time of transit of the plasma between the two stations. It 
is important, however, to note that the duration of the steady 
phase is the same whether it is. obtained from the terminal 
measurements or the pressure measurements. 

Measurements of pressure and electron density indicate 
that the one-dimensional assumption is valid as well. Pres- 
sures were measured at a station 23 cm downstream from the 
electrode insulator junction, at two representative positions 
in the channel top, the centerline of the channel and 5 cm 
to one side. The values of measured pressure fall within a 
few per cent of each other indicating a negligible pressure 
gradient across the width of the channel as expected. No 
pressure measurements were taken across the height of the 
channel. However, it is reasonable to expect the pressure 
to remain uniform across this dimension at the downstream 
position, because this position is located outside the current 
distribution at a distance which is at least three times the 
characteristic distance associated with the downstream ex- 
cursion of the current distribution. 

A typical measurement' of electron density across the 
height of the channel obtained spectroscopically is shown in 
Fig. II-3. These measurements represent an average of the elec- 
tron density along the line of sight which is the width of 
the channel. Though the measured densities show some vari- 
ation across the height of the channel, the deviations are 
found to be about 15 percent, so that for the purpose of the 
analysis electron density is treated as constant with a value 
equal to the average value. 

For the one-dimensional assumption to be valid it is also 
necessary to assume one-dimensionality in the electromagnetic 



= 22.7 cm FROM ELECTRODE 
INSULATOR JUNCTION 
= 85 kA 
= 96 g/sec 



FIGURE E-3 

AP 25- 4732 


ELECTRON DENSITY PROFILE 
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fields as well. This assumption is based upon the small 

height of the channel with respect to its width and upon 

the uniformity of the current density in the discharge as 

well as in return conductors. Previous magnetic probing 

of the discharge in the : same accelerator with similar 
70 92 

currents ' have shown that the two-dimensional effects 
associated with the finite size of the conductors are not 
too severe and that no distinct current concentrations 
or spokes exist. Therefore, a one-dimensional treatment 
of the problem is justified. 

Since viscous and heat transfer effects depend quite 
heavily upon velocity and temperature gradients in the 
flow which are not known a priori as well as on transport 
properties of the fluid, the analysis will be carried out 
on the presumption that these effects are small and can 
be neglected. It will be later shown that this is indeed 
the case. 

E. Interpretation of the Data 

The experimental data can now be interpreted within 
the framework of the steady one-dimensional conservation 
relations stated in Sec. C. The analysis will be carried 
out between two stations in the channel: 

Station 1 . - Located 13.7 cm upstream of the electrode in- 
sulator junction at the same axial position as the injectors, 
and. 

Station 2 . - Located 22.7 cm downstream of the electrode in- 
sulator junction. 

1. Momentum Conservation 

The pressure p^, density and velocity u-^ at a sta- 

tion 1 are estimated on the basis of a mass flow of 96 g/sec 
assuming choked flow at the orifices. Under these assump- 
tions: 5 / . 

jOj = 0.7 x i 0 njvn 1 - 

j> = 0.1 ^/vw 3 

U , = 3.7 x 10 Z m/s 
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For a current of 57 kA, which is the lowest current used, 

M* J Z A = 7 %\o z * » 0.35 x iO z n =f l U, Z A l 
Z w z 

7 * 10 1 n y? O.Z*iO z n = fo,A, 

9 X ID 4 /)/m z ^ 0. 5 x to* n/wi z ~ fc 

Neglecting these terms, the momentum equation becomes: 

jD 4 - f u 2 = Mo J Z ( 27 ) 

r 3 ' z** 

Table I summarizes the values of the magnetic pressure at 

station 1, / the static pressure at station 2, p, and 

zw 

the electron density at station 2, n e . All these quan- 
tities are determined experimentally. Under the assump- 
tion of full (single) ionization, the density of the fluid f 
becomes : 

f = /) e 

where m^ is the mass of the argon ions. 

There is strong experimental evidence to suggest that 
the full (single) ionization assumption is indeed correct. 
Spectroscopic observation of the flow in this accelerator 
and in a coaxial MPD accelerator operating with similar 
current densities indicate that singly ionized argon is the 
dominant species m the flow. Further, there is no 

energy source to increase the ionization level downstream. 
Extensive recombination, on the other hand, is not believed 
to take place in the bulk of the plasma since the constant 
area-channel -precludes expansion -and cooling'-of - the - flow". 

Using the measured values of p, J and n e and the sub- 
sequently calculated density f , equation (27) is used to 
determine the average velocity u necessary to conserve 
momentum across the discharge. This value is compared in 




TABLE I - MOMENTUM CONSERVATION 
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the table with the velocity obtained from a Doppler shift 
measurement u D . The agreement between the calculated 
velocity and measured velocity is good verifying the several 
assumptions involved in the calculation. 

The column rrip shows the mass flow calculated on the 
basis of the measured densities and calculated velocities 
at station 2 for the three currents studied. It is found 
that the output mass flow is a strong function of current 
despite the fact that the input mass flow is fixed at 
96g/sec. This indicates that mass is ingested from the gas 
dynamic flow established prior to breakdown, or from elec- 
trode and insulator ablation. 

2. Energy Conservation 

On the basis of the properties of the flow through the 
injectors it is found that the inlet kinetic energy is « 

negligible with respect to the input electrical energy to 
the flow, i.e. 

r.u. a «,f 

With this approximation, the energy equation becomes: 

'■fuA (f, + ft 1 ) = << VJ < 28 > 

. The enthalpy in the exhaust flow could be expressed as 
a function of density and pressure. It is, however, 1 more 
convenient to express it as a function of temperature and 
pressure for reasons which will become clear below. 

Since there are no reliable temperature measurements 
in this accelerator, the temperatures will have to be esti- 
mated on the basis of an equation of state which takes into 
account the fact that the fluid is fully (singly) ionized. 

Experiment ally, it has-been -found that- in accelerators of ■ 

this type, the electron temperature remains nearly . constant 
at 1. 5-2.0 eV over a wide range of operation conditions. 

Thus, the ion temperature can be calculated from the measured 
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pressure and electron number density, the assumed, electron 
temperature, and the following expression: 

(29) 


p = fle k T e +■ k (T e T x ) 


where T e is the electron temperature and is the ion 

temperature. Table II shows the values of ion temperature 
obtained with this equation assuming an electron temperature 
of 2 x10 4 °K. 

Consistant with the two fluid model incorporated in the 
equation of state presented above, the enthalpy can be ex- 
pressed as: 

(30) 


Here, 


T i 


$ = i ( r i, r e , f) 

describes the energy content of the random ion 


(31) 


motion while T e describes the energy content of the elec- 
tron random motion and of the excited states which in this 
case include ionization. 

The electron temperature never reaches the ion tempera- 
ture since thermal energy transferred from ions to electrons 
easily passes from the electron thermal modes into the ex- 
cited state system. This process can account for the con- 
stant electron temperatures observed. With. this model, the ex- 
pression for the enthalpy becomes: 

^ ' ' z 

The first term can be obtained directly from, a Mollier 
chart for argon while the second may be evaluated 
directly for each calculated value of ion temperature. 

Assuming an electron temperature of 2xl0 4 °K the express- 
ion for. the enthalpy becomes: 

4 - {smo 1 + S.Z2*io l Ti /.*) m*/s * < 32 > 

The column in Table II displays the value of the enthalpy 

obtained with the above expression on the basis of the cal- 
culated ion temperatures. The total power in the exhaust. 


(105) 

T e cr „ r r i 

c r dT + l jr 

'e 
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„ fuA(Vf + h u^) can then be calculated and is also shown 
in the table. This power should be compared to the total 
power available to the flow, oCV J, defined previously. 

For an MPD accelerator of comparable size and range of 
operation, the cathode and anode falls are typically of 
the order of 10 and 20 volts respectively. H3) ^ j_ s 

therefore estimated that between 10 and 20 percent of the 
input power is lost to the electrodes resulting in a value 
of oC = 0.85. The total input power available to the flow, 
0.85 VJ, is shown in the next column and is seen to com- 
pare very favorably to the total power in the exhaust. 


3. The Energy Deposition Ratio 


The initial energy deposition ratio, which is the 
central interest of this study can be calculated using two 
approaches. One approach uses the ratio of streaming to 
thermal energy in the exhaust under the assumption that the 
constant area channel does not allow acceleration of the 
flow downstream of the discharge region. The energy depo- 
sition ratio calculated by this method is shown in Table II 
in the column labeled • According to the calcula- 

tions only about 20 percent of the energy in the exhaust 
appears as directed kinetic energy. 


The other approach uses the ratio of the motional emf 
and resistive components of the terminal voltage across the 
accelerator. To define the motional emf, the energy equation 
is used in the form: 

IV - Je f Jv < 33 > 

V 

where W is the rate of energy input to the plasma by the 
external electric field and j is the current density. In 
jthe pla.sma frame, of reference,., the. electric, field -is. - - 


a' = i + u ^ a 


(34) 


where B is the local magnetic field and u is the velocity 
of the plasma. Substituting the equation above in the energy 
equation, the total energy input to the plasma by the external 
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J 


fields becomes: 


XV = $£'■ f Jv - 

V 


V 


( 35 ) 


Since the total energy deposited in the streaming plasma 
can be written as ot-VJ where as before £< is the fraction 
of the input power available to the flow, the voltage drop 
minus the fall voltages can be written as the sum of two 
voltages : 


v ' j 



- v. 


~t V. 


ud 


where V R is defined as the resistive voltage drop across 
the accelerator while V uB is the motional emf. 


To perform a rigorous calculation of V uB a detailed 
knowledge of the velocity and current density in the accel- 
eration region is required. Its magnitude, however, can be 
estimated by assuming reasonable profiles for the velocity u, 
current density j, and magnetic field B. Table III sum- 
marizes the results of calculations made for different ar- 
bitrary profiles. In this table, B-^ is the magnetic field 
upstream of the current distribution, u is the downstream 
velocity, w is the width of the accelerator, £ is the 
width of the current distribution and J is the total current. 
The table shows that the motional emf takes the form: 


V ua ~ U3,k0 (36) 

where Q is a factor which depends on the particular pro- 
files chosen and has a value that ranges between 0.16 and 
0.41 for the profiles shown. 

The ratio of motional emf to resistive drop is shown 
in the column V uB / <*rV-V uB . The value of V uB is calculated 
assuming that there is no downstream expansion of the flow. 
Therefore, the measured velocities accurately reflect the 



TABLE III - v n FOR DIFFERENT VELOCITY AND CURRENT DENSITY DISTRIBUTIONS 
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velocity of the fluid at the exit of the current distri- 
bution. The geometric factor Q for this calculation is 
set at 0.20. 

There is reasonable agreement between the energy 
deposition ratio based on kinetic energy and enthalpy and 
the ratio based on motional emf and resistive voltage drops. 
Both these ratios point to the conclusion that the discharge 
itself imparts less energy to the flow by direct body ac- 
celeration than by resistive heating. In the absence of 
downstream expansion this predicates low exhaust speeds and 
high downstream pressures and temperatures. Conversely, it 
seems reasonable to infer that the comparatively high exhaust 
speeds found in conventional MPD arcs of the same power must 
be derived in substantial part from electrothermal conversion 
from the high enthalpy arc plasma via the downstream expan- 
sion permitted in this geometry. 
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F. Appendix: Viscous Dissipation and Heat Transfer 


The calculations performed in the previous section are 
based on the assumptions that viscous drag and*heat transfer 
to the walls are not severe. The values of the flow proper- 
ties previously calculated and approximate values of the 
transport properties published in the literature can be used 
to estimate the magnitude of these effects. 


The Reynolds number of the flow based upon the length 
of the channel is, for the 85kA current pulse: 


Re 


= f u i 


where 


II-3 


= Z *‘0 


-/ 




/x 3)1 l0 ' 

/I - 2 >■ /0 


S kj /Cv\~\ 


For an incompressible turbulent boundary layer, the friction 

II-4 

coefficient is: 
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0,03 
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The ratio of viscous forces to electromagnetic forces becomes 


fv_ f fU*z£i + 
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,-2 


= 2 x /o 


s-2 


-Z 


where h = 5x10 “ m and w = 1.5x10 “ m. It is therefore 
justified to drop the viscous dissipation from the momentum 
and energy equations. 

It is not possible to perform a similar calculation for 
the heat transfer to the channel walls due to the complex 
nature of the working fluid. However, it is possible to 
estimate the magnitude of the heat transfer from the ratio 

of the power delivered to the. f.low_ to the. -total power- in- the 

f u A -HaVx.) 

exhaust which is shown in the column — — — in Table II. 

VJ 

According to this ratio, 75 to 85% of the power delivered to 
the beam appears as power in the exhaust. On this basis, it 
is estimated that 15 to 25% of the input power is lost by 
either heat transfer to the channel walls or by radiation. 
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III. QUASI-STEADY DISCHARGE PROPERTIES 

A. Time-Resolved Interferometry of the MPD 
Exhaust Plume (Bruckner) 

An extensive time-resolved interferometric survey of the 
exhaust plume of the MPD arc has been carried out to map the 
axial and radial profiles of the velocity and temperature of 
argon II, which is the spectroscopically dominant state of 
the propellant in the discharge. Detailed analysis of the 
data illuminates various aspects of the discharge, such as 
the relative roles of electromagnetic and electrothermal ac- 
celeration mechanisms, and the partition of directed kinetic 
and random thermal energies in the exhaust flow. A full 
treatment of the subject will appear in a forthcoming Ph.D. 
dissertation. 125 

Exhaust velocities have been measured in the past by 

means of displaced pairs of biased double electrostatic 

probes which detect the time of flight of local fluctuations 

85 97 118 

in the plasma number density. ' ' Unfortunately such 

probes suffer from the fact that their presence in the ex- 
haust jet disturbs the flow and they are unable to distinguish 
particle motion from wave motion. An independent measurement 
of the exhaust velocity can be made by observation of the 
Doppler shift of selected spectral lines. Although such a 
spectroscopic technique does not disturb the plasma flow, it 
is not capable of direct local measurements. Rather, the ob- 
served line shift is the integrated effect of all emitters 
along the line of sight and thus is distorted by radial and 
axial distributions of emitter number density, temperature, and 
velocity. In principle, the complete velocity and temperature 
distribution of the jet can be obtained by properly unfolding 
detailed spectroscopic observations at selected angles to the 
axis of the discharge 111 1 as long as the jet is cylindrically 
symmetric. 

In- order to achieve the high spectral resolution neces- 
sary to detect velocities of the order of 10^ m/sec, we 
have used a Fabry- Perot interferometer which can be scanned 
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electrically, via a piezoelectric mount on one of the inter- 
ferometer plates, at selected times during the arc pulse. A 
schematic of the interferometer and its supporting optics is 
shown in Fig.III-1. The output of the pulsed argon ion laser 
is superimposed on the optical train to provide a convenient 
sharp reference wavelength. The pinhole behind the projec- 
tion lens, L^, selects the center of the scanning fringe 
pattern produced by the Fabry-Perot, and acts as a spatial 
filter defining, in conjunction with the lenses ahead of it, 
the conical volume of plasma which will contribute light to 
the photomultiplier detector. The mirror M-^, can be adjusted 
by means of two micrometer screws to select for examination 
any desired line of sight through the exhaust plume. The 
entire system is enclosed in a sound absorbing box and mounted 

on a heavy table which is vibration-isolated from the floor, 
o 

The 4880 A line of argon II has been chosen for study 

because it is in a relatively clean region of the arc spectrum 

where there exist only a few relatively faint impurity 
119 

lines * and because it is readily available from the laser 
as the unshifted wavelength standard. 

Complete radial surveys have been done at three lines of 
sight which intersect the discharge axis at distances of 
9.1 cm, 15.8 cm, and 39.8 cm from the anode face, at angles 
of -38. 9°, 40.6°, and 36.3° respectively. Up to 26 radial 
positions for each axial station were examined, the MPD arc 
being fired three to eight times at each to acquire enough 
data to smooth out random shot-to-shot variations. Fig. III-2 
displays a typical data record for the 16.8 kA, 5.9 g/sec 
arc condition, at the 9.1 cm axial position, on axis. The 
spectral line is scanned several times in succession during 

the discharge, the observed fringe separation,, or .free. -.spec-. . 

- r ' " o 

tral range, being 0.74 A in this case. The observed Doppler 

shift corresponds to an axial component of velocity of near- 

4 4 

ly 10 m/sec.. Components of velocity as high as 2.0 x 10 

m/sec have been detected at the 39.8 cm location. The radial 

distribution of axial velocities appears to be relatively 
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flat, "however, a detailed map must await reduction of the 
transverse radiance and temperature profiles. 

An interesting effect observed in the velocity survey 
is the radial distribution of line radiance, which 
shows a curious dip and a second peak off-axis, Fig.m-3. 

This suggests a more complicated jet structure than has 
heretofore been assumed. 

Following the velocity measurements, the optics were 
rearranged to carry out radial temperature surveys perpen- 
dicular to the discharge axis at four axial distances down- 
stream of the anode face: 4.6 cm, 10.9 cm, 15.9 cm, and 30.6 
cm. It was during these experiments that a prominent split- 
ting of the spectral line was discovered. This is well dis- 
played in Fig. II 1-4, which was obtained at 10.9 cm from the 
anode, on axis. This effect is caused by a Doppler shift 
resulting from substantial radial components of velocity in 
the exhaust plume. This radial velocity is seen at all four 
axial stations, increasing in the downstream direction, as 
shown in Fig. III-5a and reaching a maximum of about 5400 m/sec. 
The very sharp separation of the approaching and receding 
Doppler components of the spectral line indicates that the 
emitting ions are concentrated in two distinct regions of 
the jet across its diameter. This model is supported quali- 
tatively by the radial radiance profiles seen in the velocity 
survey and also in the survey perpendicular to the discharge 
axis. A fortuitous result of this structure is that the strong 
splitting of the spectral line allows relatively unambiguous 
temperature measurement of the exhaust jet. From the Doppler 
width of the individual line components, the temperature of 
the argon II was calculated and plotted in Fig. III-5b, as a 
function of axial distance. It is seen that the ion tempera- 
ture drops quickly to a value less than 1 eV, which reflects 
efficient conversion of random thermal energy to directed 
kinetic energy of the jet. 

As an additional check on exhaust velocity, the optics 
were aimed down the Plexiglas tank to observe the discharge 



27 



o< 


30NVK3VU 3AI1V13U 


FIGURE IH-3 

AP25-4733 


RADIAL PROFILE OF LINE RADIANCE AT 4880 



28 


100 /isec/DIV— — 



THE 4880 SPECTRAL LINE OF ARGON H TAKEN 
PERPENDICULAR TO DISCHARGE AXIS 
10.9 cm FROM ANODE 


FIGURE IH-4 

AP25-P 395 


o< 



VELOCITY m/s 





30 


axially, i.e., head-on. The exhaust plume was examined at 
various radii and at several circumferential angles. It 
was found that the plume exhibits a complex, lobed struc- 
ture, with the argon II seemingly concentrated in six well- 
defined jets which are angularly in line with the six pro- 
pellant injector ports in the discharge chamber. This re- 
sult correlates well with the previous remarks on jet 
structure. The particle velocities in these jets range 

3 

from about 5x10 m/sec at the anode orifice to about 
4 

1.2x10 m/sec approximately 20 cm downstream. 

in order to verify the existence of the complex six- 

jet structure of the exhaust plume, the discharge was photo- 

o 

graphed through the 4880 A interference filter to reveal the 
spatial . distribution of .the argon II emitters. The re- 
sults, shown in Figs. III-6a and III-6b emphatically ex- 
hibit this interesting structure. The six argon II jets 
(three of which are seen in the side view) line up directly 
with the injector ports, suggesting that the injected gas 
streams remain well-collimated through the discharge. Be- 
cause this structure had not been detected in the past, and 

because spectrograms reveal an abundance of carbon and other 

105 . 

impurities in the jet, it was hypothesized that the regions 

between the argon jets are filled with ablation products from 

the Plexiglas rear wall of the discharge chamber. This is 

confirmed by Figs. III-7a and III-7b which show photographs 

o 

of the discharge taken without the 4880 A line filter, thus 
making all wavelengths visible. In the side view the jet 
structure has smeared out and is only faintly visible. In 
the end view (Fig. III-7b) it can be seen that the dark 
spaces between the bright argon jets visible in Fig. III-6b 
are now themselves luminous and extend radially outward be- 
tween- the -injectors to the chamber walls,' in a starlike 
pattern. This remarkable argon II jet structure had never 
before been detected since the various diagnostic tools 
used in the past, such as magnetic probes, electrostatic 
probes, pressure probes, etc. are unable to distinguish among 
the various ionic and atomic species present in the discharge. 



a) SIDE VIEW OF DISCHARGE AT 4880 A 
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b) HEAD-ON VIEW OF DISCHARGE AT 4880 A 


FIGURE m-6 

AP25- P 396 





a) SIDE VIEW OF DISCHARGE AT ALL WAVELENGTHS 



FIGURE m-7 

AP25P397 







33 


With the discovery of this new phenomenon the reduction 
of the integrated velocity profiles to local values must be 
considerably modified. The results of this analysis, together 
with the complete temperature and intensity profiles, are dis- 
cussed at length in Ref. 125. 
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B. Spectroscopic Studies (Hixon, von Jaskowsky) 

Continuing the previous spectroscopic studies of 

composition, temperature, electron density and velocity in 

105 109 119 

the quasi-steady MPD arc discharge, ' ' additional 

experimental results can now be presented which establish 
the presence of AIII, indicate the absence of AI , and iden- 
tify two distinct radial velocity components in the expan- 
sion of the exhaust flow. These data rely on earlier Kerr- 
cell photographs and time-resolved signals from photoelectric 
119 125 

detectors ' which show the constant level of electro- 
magnetic radiation from the discharge, thereby allowing 
time-integrated photographic recording of spectra and re- 
lieving the use of complex apparatus for very short expo- 
sures . 

a) Doubly ionized argon, AIII 

Argon III produces several strong lines in the near 

o 

ultraviolet between 3200 and 3500 A. Because these wave- 
lengths are near the absorption limits of the Plexiglas 
vacuum tank, a 7 . 5-cm diameter quartz window was installed 
in one of the tank ports. Fig. III-8 shows the complete 
arrangement of vacuum tank, quartz lens, and Hilger E-l 
Quartz Spectrograph where the data were recorded using 
Royal X Pan Sheet Film. The spectrograph, focused and 

calibrated with a mercury 8W germicidal lamp, provided a 

o 

linear dispersion between 6.4 and 8.2 A/mm over the region 
of interest. A typical spectrogram, covering an 8-cm high 
cross section of the exhaust flow centered on the axis and 
4 cm downstream of the anode, is shown in Fig. III-9. The 
spectra were compared to previous records of pinch discharges 
to identify All reference lines of known wavelengths. The 
positions of intermediate lines were measured with a com- 
parator and the corresponding wavelengths were calculated 
by linear interpolation. 

The observed AIII lines for various current levels and mass 
flows of argon are listed in Table IV. For the highest current 
level (approximately 40 kA) , at least two lines from each of 
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TABLE IV - Observed AIII Lines 


J = 40 kA m 

= 

36 g/sec 



Multiplet No. 


- o 

Measured ^A(A) 

Identified line 

(A) 

2 


. 3503.94 

3503.58 


2 


3499.67 

3499.64 


3 


3346.35 

3344.72 


3 


3336.61 

3336.13 


1 


3313.64 

3311.25 


1 


3301.89 

3301.88 


1 


3285.88 

3285.85 


J = 39 kA m 

= 

6 g/sec 



' 3 


3346.49 

3344.72 


3 


3336.64 

3336.13 


1 


3301.89 

3301.88 


1 


3285.88 

3285.85 


J = 16 kA m 

= 

6 g/sec 



3 


3358.59 

3358.49 


3 


3343.24 

3344.72 


1 


3313.72 

3311.25 


J = 12.3 kA 

m 

= 3.8 g/sec 



No AIII 

J = 8.7 kA 

ra 

= 6 g/sec 




No AIII 
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the three AIII raultiplets appear quite distinctly on the 
spectrogram. For example, two lines of multiplet No. 2 
can been seen in Fig. III-9 for the 40 kA, 36 g/sec opera- 
ting condition. At the 16 kA level, most of the AIII lines 
have faded although three lines are still faintly visible. 

At lower currents, for example 8.7 kA, no AIII lines are 
recorded and many of the less prominent All lines have also 
faded out. Although the AIII intensity obviously decreases 
with current, the total radiation from the discharge also 
decreases making it difficult to determine any change in 
the ratio of AIII to All over the test conditions. Never- 
theless, it is clear that All is always dominant for all 
data recorded. The influence of mass flow on the appearance 
of AIII is apparently secondary compared to this overall 
change in radiance with current. 

Fig. III-10 shows the radial extent of the observed 
All and AIII lines for the 16 kA x 6 g/sec discharge super- 
imposed on electron density contours derived from an earlier 

105 

study of the Stark widths of the and lines. At 

radii beyond 4 cm, appreciable vignetting occurs, and state- 
ments about the occurrence of All or AIII at larger radii 
are at best rough estimates. It is nevertheless quite evi- 
dent that the AIII lines fade out about 1.5 cm from the axis 
while All fades toward the edge of the field of view at a 
radius of about 4 cm. 

b) Neutral argon, AI 

Spectrograms of the MPD exhaust beam recorded previously 

with a glass 3-prism Steinheil spectrograph were examined 

for evidence of AI lines 1 cm downstream of the anode. Over 
o 

2000 A in the visible portion of these spectrograms, the 

position of all lines in the- vicinity of ' known AI' spectral 

locations were measured with a comparator to an accuracy of 

o 

a few tenths of 1 A. Using this technique, the absence -of 
prominent AI lines for at least six different spectral loca- 
tions has been established £br the 16 kA, 6 g/sec argon dis- 
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charge. Considering that neutral argon must exist at the 
beginning and end of the current pulse, the absence of any 
AI evidence further strengthens the experimental observa- 
tion that All is the dominant radiating state of the pro- 
pellant. 

c) Radial expansion of the MPD exhaust 

A continuing program to determine the heavy particle 
temperature in the exhaust flow by high resolution spectro- 
scopy has revealed an interesting complication. For this 
study a Fabry-Perot interferometer was used in tandem with 
the Steinheil glass spectrograph to examine a perpendicular 
cross section of the exhaust flow 10 cm from the anode. 
Figure III-ll shows a typical interf erogram of the resulting 
All lines. Many of the recorded fringes -near the centerline 

of the exhaust flow, for example the circled fringe of the 

o 

4806 A line, consist of two clearly separated fringes due 
to radial velocity components, one towards and the ( other 
away from the spectrograph. The recorded Doppler shift cor- 
responds to a radial velocity^ of approximately 3500 m/sec, 

i 

in good agreement with that radial velocity calculated from 
the piezoelectrically scanned Fabry-Perot, Doppler j shift 
measurements (see Section III-A) . Consequently, the re- 
duction of heavy particle temperatures from the Doppler 
width, as recorded by the Fabry-Perot-Steinheil combination, 
must proceed in conjunction with the reduction of the veloc- 
ity profile. A survey with a decreased free spectral range 
for even further separation of the Doppler components is 
now in progress. 
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C. Momentum Balance across the MPD Discharge (Cory) 

Measurements of the impact pressure in the accelerator 
exhaust jet, and static pressure in the arc chamber have 
been performed as part of a program to investigate the con- 
servation relations across the acceleration zone of MPD arcs. 
The pressure probe used for these measurements consists of 
a 5-cm long dielectric stub which transmits pressure, but 
not heat, electrical charge or induced EMF to a thoroughly 
shielded PZT-5 piezoelectric crystal, and a 90-cm long by 
1-cm diameter backing rod whose function is to prevent re- 
flected acoustic noise by containing the propagating pres- 
sure wave. This probe has a rise time of 10 jasec and gives 
distortion-free pressure records for approximately 1 msec. 

The probe was first installed on the moveable platform 
in the Plexiglas tank to measure the momentum flux in the 
MPD exhaust. A radial profile of the impact pressure 28 cm 
from the anode for the 16 kA, 5.9 g/sec discharge is shown 
in Fig. III-12. The profile is nearly gaussian and decays 
with the square of the axial distance from the anode. Simi- 
lar measurements for 1.9, 3.8, 23 and 36 g/sec argon over a 
range of currents from 8 to 47 kA showed this shape to be 
nearly invariant with mass flow and current. These features 
are typical of jet flow into a vacuum and suggest a model 
for the MPD exhaust once it has passed out of the acceleration 
region. 

The same probe was then used to measure the pressure in 
the arc chamber by inserting it in a suitable static pressure 
tap in the chamber side wall. Measurements were obtained 
over the same nominal mass flow and arc current range listed 
above. This measured chamber pressure was then used to calcu- 
late the electrothermal thrust component. An empirical rela- 
tion embracing 95% of the data is 

Electrothermal Thrust = (2.45 - 0.4) x 10 -5 J 1 * 5 m 0,3 


IMPACT PRESSURE , atm 
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The total momentum flux from the accelerator consists 
of both an electromagnetic term and an electrothermal term. 
For our geometry and the 16 kA, 5.9 g/sec discharge, the 
former is evaluated from the familiar self-field express- 
ion, T = u J 2 Hn (r /r ) + 3/4J/4 7 T~, to be 62 nt. The 

O ^ cl C -J 

latter component, calculated by assuming the measured cham- 
ber pressure of 0.02 atm acts over the entire anode orifice, 
is 16 nt making a total thrust of 78 nt. Integration of the 
measured impact pressure profile, which should also give the 
total thrust, yields 81 nt in good agreement with the calcu- 
lated sum of the components. Similar agreement was obtained 
at the other mass flows and currents indicating all major 
sources and sinks of momentum in the MPD exhaust have been 
identified. 
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IV. SELECTION OF OTHER PROGRAMS IN PROGRESS 

A. Hollow-Cathode, Quasi-Steady MPD Arc (Parmentier) 

It has been experimentally demonstrated that substan- 
tial benefits accrue from the use of a hollow cathode for 
currents below 10A. 2 the operating conditions 

and geometry are proper, a stable discharge originates inside 
the cathode cavity resulting in a lower cathode fall vol- 
tage, lower radiation losses, and less erosion. The object 
of the present study is two-fold: first, to investigate whether 
a hollow cathode will operate in the multi-megawatt power 
range typical of quasi-steady MPD arcs; and second, if this 
operation can be demonstrated, to determine whether the benefits 
observed at low powers are present at high powers as well. 

119 

Earlier tests with a cylindrical hollow cathode 
showed that under some test conditions, the total arc im- 
pedance decreased, while for others, the impedance increased. 

In addition, Kerr-cell photographs showed little evidence 
that the discharge was emanating from the cathode cavity. In 
recent tests, two stainless steel hollow cathodes were em- 
ployed. The first, Fig.IV-la, is a conical configuration of 
overall dimensions similar to our solid conical cathodes. The 
motivation for this design is that the arc attachment at the 
cathode tip should generate a sufficiently large magnetic 
pressure to overcome the gasdynamic pressure associated with 
the mass flow through the cathode, and force the current 
attachment back into the cavity. The second cathode, Fig.IV-lb, 
is a larger diameter open configuration with nylon covering 
both the outer cylindrical surface and the downstream face. 

With this design, the discharge is forced to attach somewhere 
on the inner cylindrical surface. 

The cathodes were tested in our standard MPD chamber con- 
figuration with the current pulse supplied by an electrolytic 
capacitor bank and the mass flow provided by the now standard 
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electromagnetic solenoid valve. A small manifold mounted 
directly behind the arc chamber split the flow between the 
six normal chamber injectors and the cathode feed tube'. 

Various orifices were inserted in these feed lines to change 
'the fraction of the total mass flow which passes through the 
cathode. , 

Voltage data for these two cathodes, acquired over a 
large range of current, mass flow, and percentage distribution 
of mass flow, showed the insulated cathode to have the greater 
effect and thus, only these data will be presented. Although 
wide in scope, the data must be considered strictly preliminary 
due to the limited number of shots taken. Consequently, inter- 
pretation of and conclusions drawn from these data will be de- 
ferred until the trends have been fully verified. 

A typical voltage-current discharge characteristic for a 
mass flow of 2.1 g/sec argon and a representative division of 
the flow of 67% through the chamber injectors and 33% through 
the inside of the cathode is shown in Fig.lV-2. This type of 
characteristic has been observed over the complete spectrum of 
flow division percentages for the same total mass flow, a sum- 
mary of which is shown in Fig.IV-3. Also shown in the figure is 
one possible set of lines through these data for fixed current 
levels for the purpose of representing performance trends. For 
a fixed terminal voltage and a fixed flow division between 
chamber injectors and cathode, the measured current as a function 
of argon mass flow is shown in Fig.iV-4. The appearance of a 
maximum current followed by a plateau as mass flow is increased 
is borne out for other voltages and flow divisions. 

These trends are summarized in Fig. IV-5, where the mass flow 
is plotted against percent division of the flow for several 
values of current and a fixed terminal voltage of 220 volts. 

For a fixed percentage flow inside, the data indicate that there 
are two different injected flow rates where the power input into 
the discharge is the same. Stated in another way, there is an 
optimum mass flow rate (defined as the maximum current for a 
fixed voltage) for any given distribution of the flow between 
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outer injectors and the cathode. The locus of these optimum 
mass flows coincides with a line of constant J 2 /m, the con- 
stant equal to 75 for these particular data. For comparison, 
a line of J /m = 52 has also been drawn which is the charac- 
teristic value calculated for this geometry from a simple 
argument based on an equipartition of energy between kinetic 
and thermal modes. a similar line drawn for J 2 /m = 52, 

where m is only that mass passing through the cathode, seems 
to bear little relation to any of the key features of the 
figure emphasizing that the total mass flow is a better- scal- 
ing parameter for these terminal data. 

Data points which lie above the J 2 /m = 75, i.e. mass flows 
greater than the optimum value, may indicate an overfed region 
whereas those below the line may represent a region of starva- 
tion. Thus, the apparent insensitivity of the data to the 
initial flow distribution at low mass flows may be a manifes- 
tation of electrode or insulator erosion at these conditions. 

In support of a starvation argument, the front face of the 
cathode and insulator were severely eroded after this series 
of shots. 

To minimize this erosion and to investigate these trends 
in more detail, a new tungsten hollow cathode has been in- 
stalled in the discharge chamber. The cathode is larger, 3.2 cm 
in diameter, to facilitate probing of the interior in order to 
determine what fraction of the discharge current is emanating 
from the cavity, and to allow better photographic observation. 

In addition, this configuration will be driven by one of our 
high voltage capacitor lines instead of the electrolytic line 
to insure a flat current pulse at a predetermined level. 
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B. Time-Resolved Thrust Measurement (Saber) 

Recent experiments have shown the first positive results 
of a thrust-time profile measured with piezoresistive strain 
gauges mounted on structural elements of a quasi-steady MPD 
accelerator. Initial tests revealed that the response of 
individual gauges mounted on the inside of the cathode cylin- 
der was dominated by both electrical and mechanical noise. 

The electrical noise is now reduced to a tolerable level by 
the various techniques illustrated in Fig. IV-6. Most of the 
electromagnetic pickup is eliminated by differentially moni- 
toring the output of two identical strain gauges, one mounted 
on the cathode cylinder and a second passive gauge suspended 
directly above it. Careful shielding of this sandwich struc- 
ture as well as the rest of the bridge circuit leads to a 
noise level of 0.1 mv compared to an anticipated signal level 
of the order of 1 mv. Mechanical pickup by the sensor was 
primarily related to initiation of the mass flow in the shock 
tube mass injection system. This source of noise was tem- 
porarily eliminated by operating in the "ambient" mode, i.e. 
prefilling the entire Plexiglas vacuum tank to a fixed pres- 
sure level, say 100 p, and firing the accelerator without fur- 
ther injection of any mass. This problem will be permanently 
resolved when the electromagnetic valve mass injection system 
is in standard use. 

Fig. IV- 7 shows the response of the strain gauge to a 40 kA 

x 0.5 msec current pulse on three different time scales. The 
50 psec delay before the initial gauge response is due to the 
finite propagation time of the strain wave down the cathode cylin- 
der, but also illustrates the low noise level compared to the data 
signal. After 50 psec, the trace displays a superposition of two 
characteristic frequencies on the data signal - a 10 kHz pattern 
corresponding to the propagation and reflection of the longi- 
tudinal strain wave in the cathode cylinder, and a 1 kHz wave 
corresponding to a "drumhead" oscillation of the cathode end 
plate. These oscillations are more clearly shown in Fig. IV-8 
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where the gauge response to both a forward and reversed 
21.5 kA x 1 msec current pulse is shown. Since the thrust 
is almost completely generated by the interaction of the 
current with its own magnetic field, current reversal should 
have no effect on this electromagnetic thrust component as 
long as the current distribution remains azimuthally sym- 
metric. As seen in Fig. IV-8 there is no discernible dif- 
ference in the transducer output for these conditions indica- 
ting that the observed signal is both free of spurious noise 
and representative of an electromagnetic thrust record. 

Further evidence that the gauge is responding to electro- 
magnetic thrust is presented in Fig. IV-9. As the current is 
raised by approximately 75%, the amplitude of the initially 
excited strain waves increases by approximately the square 
of this current change. Thus, although the mechanics of the 
accelerator deformation complicate the extraction of the 
true thrust level, the strain gauge technique appears to be 
a viable one for recording thrust-time profiles in this hos- 
tile environment. 
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C. Time-of-Flight Velocity Probes (Boyle) 

Signals from double electric probes biased to draw ion 
saturation current have been previously used in this labora- 
tory to determine plasma velocity profiles in the exhaust of 

85 97 118 

an MPD accelerator. ' ' If one assumes that the micro- 

structure of the velocity probe signal is characteristic of 
fluctuations in local plasma properties and that these 
fluctuations convect with the streaming motion of the exhaust 
flow, then the plasma velocity can be readily obtained from a 
time-of-f light . measurement based on correlation of the sig- 
nals from two of these probes displaced a known axial distance. 

A program has been initiated to investigate the accuracy 
of the calculated plasma velocity based upon the time- 
of-f light hypothesis. Theoretical support for the program 
is based on collisionless probe theory for the cylindrical 

IV-4 

Langmuir probe in a stationary Maxwellian plasma by Laframboise. 
When such a cylindrical probe is aligned at an angle to the in- 
cident plasma flow, the probe current is found to be a function 
of local plasma properties, plasma velocity, probe bias, and 
ratio of probe radius to Debye length. Since the local elec- 
tron number density, electron temperature and ion temperature 
have been previously determined in the MPD exhaust , ^ ^ 

such a probe should yield the plasma velocity from the experi- 
mentally determined voltage-current characteristic and probe 
geometry. Thus, a time-of-f light probe consisting of two of 
these cylindrical probes provides data for three independently 
calculated velocities - two from the individual probe charac- 
teristics via probe theory and one directly from convective 
fluctuations. Agreement among these velocities establishes 
the validity of the time-of-f light hypothesis on a physical 
basis . 

As a first step in this program, a double probe has been 
constructed from 0.076-mm diameter tungsten wire. The elec- 
trodes are each 7.6 mm long with a 1 mm separation. For reason- 
able estimates of local plasma properties, the electron-electron. 
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electron-ion and ion-ion mean free paths are greater than 
the radius of the probe electrodes implying a collisionless 
flow regime. A length-to-diameter ratio of 100 assures end 
effects are negligible . IV 5 The dimensions of the previously used 
parallel plate boundary layer probes were such that the probe 
operated in a transition flow regime, a difficult one to de- 
scribe analytically. 

Initial tests with this probe in the MPD exhaust indicate 
a similarity between the recorded currents of the collision- 
less probe and the previous parallel plate probe. For example, 
arrival times for the initial plasma front of approximately 
8000 m/s are recorded, and. the later quasi-steady portion of 
the signal has high frequency current fluctuations superim- 
posed on it enabling convective fluctuation measurements to be 
made. As a result, a second, cylindrical double probe has been 
constructed and voltage-current characteristics are now being 
recorded to evaluate the plasma velocity. 
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D. Repetitive Pulse Propellant Injection (White) 

As the quasi-steady propulsion concept continues to 
evolve toward a working space thruster, increased attention 
must be directed toward various other components of the com- 
plete thruster system. The two areas of most concern in- 
volve the power preparation and gas injection subsystems. 

Of these, the conversion of DC solar power into intermittent, 
high current, square pulses is reasonably straightforward, 
with the primary task that of improving capacitor quality. 

On the other hand, the conversion of a high pressure , stored 

7 8 

working fluid into 10 - 10 discrete mass pulses with strin- 

gent requirements on pulse rise time and shape, and minimum 
leakage between pulses has less evident technological solu- 
tion. Conventional valves are too slow; valves that are 
overdriven to improve their rise time have questionable life- 
time and reliability. In an effort to initiate some interest 
in this aspect of the problem, and to contribute to its solu- 
tion, we have undertaken a program of repetitive-pulse gas 
injection which explicitly eschews reciprocating valves. 

The most promising concept identified so far replaces 
the impulsively started and stopped member in conventional 
valves with a purely rotational system running at constant 
speed, where the seals will absorb far less energy, hence 
wear less. A prototype of this turbo-gas injector scheme 
has been designed, constructed, and tested which can deliver 
argon pulses of several grams per second, in a range of pulse 
lengths and separations. The device, shown in Fig. IV-10 con- 
sists of two coaxial slotted discs housed in an appropriate 
support and manifold assembly, and rotating at different 
angular velocities. Argon propellant, fed into a small 
plenum directly behind the discs, passes through the valve 
only when the slots become coincident. The gas pulse density 
profile in time is determined by the absolute speed of the 
fastest disc A in Fig. IV-10, which is coupled directly to a 
conventional electric motor; the dead time between gas pulses 
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is fixed by the relative speeds of the two discs, with the 
slower disc B geared down from the high speed shaft. Thus, 
varying the motor speed, gear ratio between the discs, and 
number of holes in the slower disc allows this system to 
satisfy a wide range of pulse demands. For example, the 
lower line in Fig. IV-11 shows the anticipated operating range 
of pulse length and pulse frequency available to the proto- 
type valve for a fixed gear ratio of 24:1 and a single slot 

in the slow disc. Operation along the line is possible by 

-3 

varying the speed of the driving motor. The 10 duty cycle 
(ratio of on-time to total on-time plus off-time) , which 
exists for all operating conditions along this line, is a 
consequence of the gear ratio, slot dimensions, and the single 
slot in the slow disc. Other operating lines accessible by 
increasing the number of slots in the low-speed disc are 
shown in the figure along with their respective duty cycles. 

For the initial testing of this valve, the rotation rate 
was set at 1725 rpm which should ideally produce 0.83 msec 
pulses with a risetime of 83 ;asec, at a repetition rate of 1.2 
pulses per second. Figure IV-12 shows the density profile for 
these pulses as recorded by a Millitorr fast ionization gauge 
a few cm in front of the exit orifice. The lower trace shows 
that a single pulse occurs approximately every 0.8 sec as 
expected with no significant leakage during the dead time. 

The detailed profiles of each of the six pulses shown in the 
lower trace are overlaid in the upper trace on an expanded 
time scale. Although the starting transient is not ideal, 
the bulk of the pulse has a risetime of roughly 200 ^isec, 

These reproducible data are highly encouraging in that they 
indicate rise and fall times and pulse durations close to 
design values and considerably shorter than hybrid solenoid 
valve systems. 

The ripple or "false start" observed at the beginning 
of the pulse is apparently caused by improper initial align- 
ment of the two discs. To improve this overall starting 
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transient, a mechanical linkage will be added to the valve 
to allow continuous adjustment of the disc coincidence dur- 
ing operation. In addition, both discs will be Teflon coated 
to reduce friction and wear. Further testing should indicate 
the full range of deliverable gas pulses, long-term wear rates, 
and the compatability of this system with the discharge. 
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B. 
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C. 

Technicians 

D. 

Supporting Staff 

Employee Benefits 


(19% of IA, IC, and ID) 
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Total Direct Costs 
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5,764 

0 
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Indirect Costs 
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TOTAL 
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